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-ABSTRACT 


Theory is presented on the lift and pressure 
distribution of bodies, wings and wing-body combinations at 
subsonic speeds. It is found that ■varying source ring 
distribution over body helps to simulate the wing-body 
interference effects. To study the phenomenon of wing-body 
interference, various mathematical models representing the 
given configurations are studied. 

The design of an aero dynani c al ly compensated 
pitot-static tube is made to nullify the position error. 
Aerodynamic compensation and insensitivity to angle of attach 
are che important features of designed pitot-static tubes. 
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PART - A 


STUDY OE WING BODY INTERFERENCE 



CHAPTER - 1 


INTRODUCTION 

i i 1 The presence of wing around body* puts the body 

at an angle of attack to an otherwise undisturbed stream of 

fluid parallel to the tody axis. This change of direction 

of fluid around body results in the wing interference* To 

simulate this effect the body can be ■ represent ed by a varying 

source ring distribution. The wing-body interference effects 

become more and more predominant with increasing yawing and 

pitching moments and angle of attack. The flow around a wing- 

body combination is studied for zero angle of attack. Por 

other angles of attack the body starts shedding vortices of 

comparable strength over the wing vortices. The modified 

wing vortex distribution changes entirely when the 

presence of body vortices are considered. The flow being at 

a zero angle of attack, the effect of vortex distribution on 

1 

body is considered negligible (Nielsen ). This assumption 
is not valid for high angles of attack. The effect of body 
vortices must be taken into account for a>10°. The effect of 
body shedding vortices has been compared with the low angle 
theory and the verification has been made by the experimental 
work carried out by H.Schliehting 2 . To simulate the effect 
of wing interference, we rely ' upon varying source distribution 
since the angle of attack is taken to be sere. 

A study of wing-body interference has been made 
for several wing-body configurations. 



CHAPTER - 2 


MATHSSIAT ICAL P ORMULAT ION 

2 j i in order to study the flow field of wing-body 

combination, the wing* body and wing-body Combination is 
represented by a set of singularities* The strengths of 
which are found by satisfying the tangential flow boundary 
conditions at the selected control points to get the given 
Configuration, The shrouded part of the wing is Considered 
•to bo of Kero camber and zero thickness* The theory applies 
to both subsonic and supersonic flows, since the linearised 
potential function and Ptardtl-Glauert equation hold for these 
flews. The similarity solution of Yon Karma n may be used 
to extend this theory for transonic flow. 

2.2 Wi ng a lone t 

The potential function for the wing satisfies the 

linearised Prandt 1-Glau ert equation^ - , 

W w W 

(1-Iv! 2 ) $ + +( J ) =0 2.1 

The failure of lifting line theory for swept ba.ck 
wings lights up the idea, of surface discontinuties . . The wing 
plane is assumed to be a horizontal plane. This plane is 
divided into a large number of panels. Over each panel a 
singularity of constant strength is distributed. The 
dcwnwash at each control point is determined. The strength of 
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eacii panel pressure discontinuity is determined to satisfy 
the flow tangency boundary condition required to obtain the 
given camber line. 

!Ehe expression for potential function for constant 
pressure singularities is given by^ 


$ 


W 


E z(x-^ ) 

* //( 


£ ) r(x-j7)io-'^Vr(y- n ) 2 +z 2 ; 


dt v 


(1-E) 


% 


ss 


z dX d'A_ 
T 1 

V a. ) +z 


- 2.2 


where E = 0.5 for M < 1 
and E = 1 .0 for M > 1 


x 



y = y p - y k 


z = 



The limits of integration are defined by the panel 
boundaries, For one corner element the potential function 
is given by 

. W _ K a u _ 0 o 

™ k ~ 2'i"~ 1 ( x-ly ) . (F3+P4 )+z^(1+1-M^). 

F2-1(F1~F5 )— yF6 v | 

f j 

where 2, p v u = pressure discontinuity across panel 


2.3 
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56 = ?+d for m<C 1 

r 

= for M> 1 


r = f~ 

y 2 + 

~' a 2~ 

d = 

-M'2y r 
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+ (1 

9 v 

x = Ik 

-M'" )y 
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- I>y 
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r “ ^(x - Iy) 2 + {-£■ + 1 -M 2 ) z 2 


d = f x «2 + r «2 


by 


The potential at a point P (x^,y‘ , z ) is given 


y v W W 7/ W 

$" (?) = * r 4> + <t> 


2.4 


where (!) W K = 1 4 

K J 

is the potential due to one corner element of constant 
pressure discontinuity: 

The perturbation velocity components are given by 

2.5 


a ^ 

U = fcf- 


V 


- 14 

= By 


9 (j) 

w = aT 


The partial derivative of eqn. (2.3) gives the perturbation * 
velocity component due to each corner clement of constant 
pressure panel. 


K '■ u 



2 TC 

K A33-, 

2 7t 

J.U.. 

2 7t 


( F3 + P4) 

1 (P3+P4 ) -z.F6 

!~(l 2 + 1 - M 2 ) P2 - 1 (P1-P5)-y.P6 


] 

) 

) 

) 

) 

) 


. 2.6 


The perturbation velocity components due to vortex 
panel at any point p(Xp,y^j z p ) are given by 
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u(p) = u 1 -u 2 -u 3 +u 4 

v(p) = v i" v 2~ v 3 +v 4 
w (p ) = w r w 2 *-w 3+ w 4 


w(p) = slope at the control points (p) x free stream velocitjr. 
This gives a set of n simultaneous linear equations for n 
control points. The solution of these equations gives the 
vrhLue of pressure diseont inunties across the panels. 



c = c - c 

P -P lower ^ upper 

4u 

or 0 p =2 .(— TT~) 2.9 


To study the pressure distribution and lift 
distribution of the wing, three configurations are chosen. 
A tapered wing, delta wing and constant chord wing have 
been treated with the described mathematical formulation. 
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Ike lift and pressure distributions are plotted and compared 
in figs. (1-15 ). Ike details of tke wings dimentions are 
given below: 


WING 

A 

s 

c 

r 

°t 

T.R. 


Tapered wing 

2.9 

10 

5 

2 

2.5 

45° 

Delta wing 

2.9 

10 

5 

0 

CO 

-d o- 

Constant chord wing 

4 

10 

3.2 

3.2 

1 

-do- 


Tke pressure distribution increases toward tke tip 
of tke constant chord wing. This is because of tke leading 
edge sweep. Tke center of pressure, therefore, moves away 
from tke root section. Tke delta, wing has tke nearest center 
of pressure to tke wing root, as compared to other wings. 


2 .3 Body, alone: 


Tke potential function for body satisfies tke 

4 

linearised Prandtl-Glauert equation 

B B 


(1-M 2 ) 




B 


4* 




(j) + 

77 


0 


. . . 2.10 


zz 


Tke body is. represented, by a constant source 

sink distribution. A source ring is a distribution of 

5 

sources of constant strength along a circle. 

The velocity induced by clement of strength q.ds 
of a source ring of radius r at a point (x,r,f> ) is given 
by fig. (1-f) 

q. ds R 
dv = r> — ~ 

4m rR r 



Taking the integration around the circumference of the ring 
we obtain for the velocity components in cylindrical 
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coordinates . 


v fe r) = 1 


2 % 


X 


9 '■*’*< " 4*r’ i 


and v^Cx.r) = ~r“ £ 


bhrhl-2r Cos (4~4 1 ) 
1 r-Sin d> 1 


— — — .. — — . — ~~ 3 

' x 2 +r^+1-2r Cos (<H l ) 


.. . 2.11 


where x = x/r* 


r = r/r ( 


for simplicity only the x and r are replaced by non-dinensionalized 
values of x and r 

r*= source - ring radius 

x = distance of control point from source ring 
r = distance of control point from x axis. 

These integrals are reduced to complete elliptic 
integrals. The expression for velocity components are given 
by the followings s 


2x 


2td? * r y ( r+ i ) z " 1 %*+ (r-1 ) 


E(k) 

I fm 


V < x - r > = 
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v 


(1- ... . ,)E(k) j 

x Z +(r-1) 2 _■ 


wliere K(k) = £ 


n/2 


lTll?sIn 2 a 


da 


E 


.V 2 p p 

(k) = ^ V1-k^ Sin 2 a da 


-D(k) = £ 


p/2 Sin 2 a 


V~ l-k 2 ''sin 2 ~a d<X 


K-E 

i-2 


and k = 


4 r 

'o' “ » — y 

x + ( r+1 ) 


2,12 


2.13 


To simulate the effect of body volume, eighteen 
source rings arc chosen with the control points placed between 
successive sources. The boundary conditions are satisfied 
to obtain tangential flow over the control points . Thus 


v ( x > r ) 

-H s = slope at any point p (x,r; 


Three bodies, an elipsoid of revolution, a. 
Sears-Haack body and a Karman Ogive of the sane length and 
volume are represented by this type of mathematical model. 
Their pressure distribution and the strength of equivalent 
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source-sink distribution is compared in figs.( 33 ). . 


The contours of the three types of bodies are 
expressed by the following equations. 


1 . Ellipsoid of revolution 

2 2 

= ” 1*0 with a/b = 10/1 

a b 2 

2. Sears-Haack body 

S (x ) = J„ 6 ; Volume (l _ J*-)2 " p /2 

3id ^ 1 


.2.14 


...2.15 


where S(x) = X-sectional area of the body at any distance 
x from the origin 
1 = length of it 's axis 

4 . 2 

volume = ^ a 7 t b ( volume of ellipsoid of revolution) 


3 . Ear nan Ogive 

)= . S i 1 2 (71-6+ 2 Sin2e } 
7C 


s(0/ 

where 

S(l) = maximum X-sectional area 


2.16 


0 = Oos Vi- , 

1 
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2 . 4 Wing-Bo d y_ Comb inations 

For the combined flow field of a wing body 
combination, the total velocity potential is not a simple 
sun of the potentials of the body <J) B and the wing $ W alone. 
The resultant flow should satisfy the boundary conditions on 
both wing and body together. Bor a wing body combination, 

the resultant potential function may be written as 

B W i 

$ + ^ + $ ..2 

Many approximate mathematical theories have been put forward 

for the explicit expression of the interference velocity 

potential function : ’ Interference potential can be 

expressed as algebric sum of interference potential due to 

1 2 

wing on body and due to body on wing . So the expression 
for (j) 1 can be written as 

(j) 1 = <j ) Wi + $ Bi 

■ 0 

The expression for T M through this equation is written as 

t M = iJ> B + <f + 1 + 4> B1 

/ + $ B1 ) + ( 4 B + $ Wi > 
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where 




]0 


and 


$ 


= $ ff + <) B1 
m = + 


2.18 


of body a 


is the potential function of wing in the presence 
3 $ J is the potential function of body in the 


pres encc 


£ wing. These potential functions satisfy the 


linearised 

except ^ 


prandtl-G-lauert equation for all Mach numbers 


(i-m 2 ) <£ + <»* + 4 


m 

r zz 


= 0 


2.19 


o -*) c + c + = 0 


2.20 


function 


tpo determine the expression for velocity potential 
ih© wing is replaced by the pressure discontinuities 

» 

horizontal panel and the body is replaced by a 


across th® 

1* tributi 011 source ring singularities. 

rphe potential function for constant pressure 


singul a ri 

component 


- t fcs is given by eqn. (2.2) and the vertical 


0 f perturbed velocity by eqn. (2.6). The potential 


for the 


and the induced velocity may be written in the 

m - m M 


fo® ^(p) = -:*T ' V.t 4 * 



The downwash at any point P( x , y ,z ) is given by 

P P P 

w(p) = w 1 - w 2 - 

= C A7* u n ^' . 2,22 

where C ¥ = IJ(l< 2 + 1-M 2 ) , 12-1 (11-15 ) -126 ~' k=1 

- | (L 2 + 1-M 2 ) .12-1(11-15 )-yl6 l 

L - Jk=2 

- <:l 2 + 1-M 2 ). 12-1 (11-15 )-126 ) 

L J k = 3 

+ .- ! l 2 +1-M 2 ). 12-1 (11-15) -H6 l ~j 

i. Jk = 4 — ...2.23 

No explicit analytical expression can be given for the 
velocity potential function of a varying source ring 
distribution, however, numerical values can be obtained by 
graphical integration of velocity field. 

When the axis of body of revolution is inclined to 
free stream the varying source - ring distribution is taken 
as singularity. This simulates the effect of wing 
interference on body. The use of additional vortex 
distribution on body is not required, since the problem is 
studied for zero angle of attack. Only at high angles of 
attack, the body starts shedding vortices of strengths 



14 


comparable to the wing vortices. This modified theory for 

high angle of attack is discussed Toy IT els on”* . The potential 

theory witl; vortex distribution on body is in good agreement 

with low angle of attack theory. The departure of the two 

curves for normal lead distribution starts at an angle of 

attack of about 10°. Ref. 1 . 

The expression for velocity components of a source 

5 

distribution of varying strength on a circle are given 


below: 

“x (x 


q. m 2x (2K(k)-2I> (k)-E(k) ) 

»r»<*) = *: Sin< i> 


2TCC ’ (x 2 +(r-l) 2 ) fx 2 +(r+l) 2 

Sin2 4 


V (x,r,4 ) = 
y 27ir (x^+ (r-1 ) hfx 2 +( r+1 ) 2 


! ( r ~ 


-)(2K(k)-2D(k)-E(k))+8(E(k)~D(k))~E(k)) 


w 2 (x,r,4) = 


^MB . 

2itr 1 (x _2 + (r-1 ) 2 ) fx 2 + (r+1 ) 2 


j ,2r Sin 2 ^ + Sin 2 cj>)’ (2K(k)-2D(k)-E(k)) 

-8 (1-2 Sin 2 4) (K(k)-D(k) )-2 Sin 2 (| E(k) j 


2.24 
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and 


where 


w z ( x >r,4 ) = Db 9^ 

Sin£c|> 1 

‘ SST S^iWy^T) 2 


iW S 2 ){(2K(k)-2D(k)-E(k))+8(E(k)-D(k))-E(li;)] 
I ^ k J 


2.25 


The resultant vertical conponent of velocity due 
source-sink ring and vortex distribution is given by algebric 
sun of the two velocity components. 


For n pressure panel singularities find n source ring- 


singularities, there are n control points on the wing and 


n control points on the body. The total down?/ ash at any control 
point of the wing is given by algebraic sun of the velocity 
due to wing pressure singularities eqn. (2.22) and those due 
to source ' ring singularities (eqn. 2. 25) 

Y\ 

w ?/ = ^ iA%W s i + Jh % j to, 3 .2.26 

j=. s d 

The dov/n?/a,sh at any of the body control points is given 


by 


n 

Wt, = 2 


B - Jl %B,i £ %w>1 1 ^ . 


....2.27 



The downwash on (n+n) control points is given by 




n 

■. z > 
1 = 1 


,i + ^ 1 B B, j 

cJ * 


2.28 


w 


-S-— = surface slope at the (n+n)- control point. 


The system of simultaneous linear equations is 
solved by Gauss -Jordon method to satisfy the tangential flow 
conditions. The final values of constant pressure 
discontinuties and varying source ring distribution are 
obtained. The new potential function for the wing in the 


presence of the body is given by 

. _ K&Urjf r 

V ' 1-: “ ' (x-ly).(F5+P4) 

+ z \ (I 2 +1-M 2 ) T'2-l(P1-S l 5 )~yP6 "Y ! 


2,29 


where 2 ,.\u is the pressure discontinuity across the panel. 
The value of new potential function at any point 


2( x p » y p? Zp) is given by 
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The value of the interference potential due to body on wing 
is given by cqn/2.18, thus at any point P 

<} 31 (P) = $ m (p) - f ( P ) £<: 


The expression for the body potential due to the source ring 
distribution is given by 






v „ rex: 
qr 


2.32 


U JL. C-u- 


The expression for peter. 


ring distribution is given by 


^ u ^ + j ^ y rly + / w 


o 


2.33 


The expressions for v v . u, v, w, are 'given by 

qjL qr 

B MB 

eqns.(2,12 and 2.24) so the expressions 'for ($ and ($ 


are known. The interference potential due to a wing on a 
body is defined by eqn 42 . 18 .) 


f 1 = f 3 . <), B 

x y z 

- S u dx + / v dy + / w dz 

o o 0 


X 


/ V r dx 

o 


r 

/ v dr 
o qx 

. . . . .2.34 
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The mathematical model described to simulate the 
given configurations can work upto + 10° of angle of attack* 
Beyond this range of angle of attack the shed< vertices of 
body are to be taken into account. To take care of body 
vortices, the body in the interference region is divided 
into a large number of panels. Each panel represents a 
constant pressure singularity, The tangential flow 
boundary conditions are applied to determine the strength 
of these singularities. 



CHAPTER ~ 3 


RESULTS AND DISCUSSION 

The study of wind-body interference involves the 
importance of satisfying the tangential flow boundary conditions 
both on wing and body. The effect of interference due to body 
on wing is high in the vicinity of body only. The wing 
interference is found to be predominent over a long distance 
on the body. As the length of the nose of the wing-body 
configuration increases, the center of pressure moves 
outward on the wing, therefore, a short nose-body should be 
preferred for high nanoeuverability . The combination of 
elliptic body, Sears-Haack body and Karman Ogive body with 
tapered wing, delta wing and constant chord wing has been 
studied for three locations of wing on the major axis of body. 

Eor a long nose shaped constant chord wing-body 
combination, the center of pressure lies outward to the 
point of center of pressure for a delta wing-body combination. 
The change in lift distribution of the constant chord wing 
due to the presence of body is not effected much as compared 
to the change in lift distribution of tapered and delta wings. 
The change in source sink singularities and wing pressure 
distribution is found maximum for a delta wing and Sears Haack 
body combination. The change in lift distribution of wing due 
the presence of Karman- Ogive body is least as compared to the 
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presence of elliptic or Sears -Haack body. The magnitude of 
interference pressure coefficient due to a delta wing and 
Sears-Haack body combination is maximum than the value of 
interference pressure coefficient due to any other combination. 

Short nose elliptic body and tapered wing combination 
is found to be the best case for least interference. The 
center of pressure of the wing, therefore, does not shift 
much due to the presence of the body. 

The wing-body interference phenomenon can be 
studied and the values of the interference pressure coefficients 
and lift distributions can be found through the computer 
programming given in appendix A. The values of pressure 
coefficient of a wing-body flow field may be found at any 
location. The value of pressure coefficient may be used for 
the design of an aerodynamic ally compensated pitot-static tube 
described in Part B, 



PARI - B 


DESIGN OF AN AERODYB AMICABLY 
COMPENSATED PITOT -STATIC TUBE 



CHAPTER - 4 
INTRODUCTION 

4*1 Pitot pressure and free stream static pressure are required 

to determine the altitude, rate of climb and air speed &f the airplane. 
The presence of a static port in the disturbed flow field created 
by the airplane causes it's port to sense inaccurate static pressure* 

To meet the need of the present day busy airports and high 
performance airplanes it is very important to know the accurate 
speed and altitude. The error in static pressure sensing is a 
'function of Mach number, angle of attack and change of configuration 
of the airplane. The' change in configuration of the airplane 
means the operation of flaps, elevators or ailerons. The operation 
of flaps can cause a non-negligible change of static pressure at 
the static port location. The error in static pressure sensing 
port due to the perturbed field of the airplane is known as the 
position error. The design of the pitot-static tube is so made as 
to provide a negative compensation to the flow field perturbation 
created by the airplane. This is known as aerodynamic compensation 
for the static pressure sensing port. The data for the design of 
the aerodynamic compensator is obtained from Part A. The 
aerodynamic compensator is designed to sense the accurate static 
pressure over a large range of Mach number. 

The design of a pitot static tube for sensing accurate 
static pressure must take into account the static pressure variation 
with Mach number and angle of attack. These static pressure ; 
variations can be minimised by a suitable location of the static 
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pressure sensing port in the vicinity of the airplane. The best 

location for the pitot-static tube mounting is found to be the . 

1 

nose boom, ahead ox the body of the ainpiane, since the variations 
of the static pressure perturbations with Mach number and angle of 
attack is a minimum at this location. The nose boom location for 
the pitot static tube has disadvantages of ma intenance and storage 
and the increased length of the aircraft complicates the ground 
movement problems . 

Since the effect of the upwash field of the wing is 
predominant in the vicinity of the wing, therefore, the second 
preferred location of the static port f s found to be the fuselage 
dorsal mounting. The third preferred location on high performance 
airplanes for static pressure sensing port is the wing tip. 

A pitot-static tube design has been tailor made to 

sense almost the exact free stream static pressure over a large 

range of Mach number and angle of attack. The designed contours 

of the pitot static tubes, shown in figs.( 55,57, .59) are 

aero dynamically compensated for the Mach numbers ranging between 

(0.3 to 0.85). The linearised theory used here can not be 

applied for MC0.3 due to the high angle of attack involved. 

Por 0.3<M<0.9, the error A,C in the cancellation of the static 

p c 

pressure perturbation is less than O.075 of C the 

Pc 

compensating static pressure. In the following Chapters of 
Part B the procedure for the design of aerodynamically 
compensated pitot-static head is discussed for the three pitot 
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head locations. For this purpose an elliptic fuselage with 
a 45 swept wing centerly mounted on it is assumed. Other 
wings and fuselage shapes and positions may be considered on a 
same procedure to design the pitot*-static tube. 



CHAPTER - 5 


DESIGH POR AERODYNAMIC COMPENSATION 

5.1 Description of ^ the_ aero dynamic ^ compensator? 

In Part A the static pressure perturbations;, 
were calculated for several combinations of wing -body shapes 
and locations. The perturbation field created by the fuselage 
and the wing, is plotted against Mach number for the three 
positions of the pitot-static tube viz.(i) ahead of body 
(ii) dorsal side of the fuselage (iii) ahead of the wing tip. 

The contour of the pitot-static tube is now sought to develops a 
negative static pressure perturbation at the location of the 
static pressure sensing port to canel the static pressure 
perturbation due to the wing-fuselage combination. The static 
pressure sensing port is made insensitive to angle of attack. 

Dor this purpose the static pressure sensing ports are located 
diagonaly opposite in the lateral direction. The shape of the 
aerodynamic compensator is made smooth to avoid abrupt 
compression and expansion of flow. We can also determine the 
perturbation field created by the airplane at any location using 
numerically the computer programme given in appendix A. 

Appendix B gives the computer programme for obtaining the 
nose contour of an aerodynanically compensated pitot-static tu.be. 
Giving the length of the pitot-static nose contour and pressure 
coefficient required for negative compensation as the input da.ta, 
the desired compensating shape of the nose contour can be obta.ined 
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We have the Bernoulli's equation for compressible 
steady flow between the free stream and the point of location 
of static pressure sensing port. 

• Po , _£L_. 5.1 


i^ + 


K + 


The local static pressure is expressed as a function 
of the free stream static pressure by 

■ p = P 0 '-(i+e),. 5.2 

' -s * 

, • Bor an' undisturbed flow the perturbation velocities 
u, v, w are zero and hence 6 is zero. Thus in the undisturbed 
flow we can sense the correct static pressure. 

. Combining ..eqns:." (5 .1 ) and (5.2) we have 

p. " i .2 v p o(i+e) 


I- ^ 


y \ 

~ 


12 

0 = 2 7 + 


.... .5 .3 


where V T (u+u)^ + -y 2 + w^ 


from which 6 can be written as 


e = - 


• 1 2 
“ L uu + ' 2 ((1-M 2 ) u 2 +v + w 2 ) I 


,5.4 


Another expression for 6 fallows from eqn. (5.2) as 


P - P, 


Defining the pressure coefficient as 

P - P, 

2 1 


°P = 


o 

tT 


.5.5 



27 


The' expression for pressure coefficient can be written as 

-[Uu + jr ((1-M 2 ) u 2 +v 2 +w 2 )~j 

— — — v. • , 5.6 

Combining oqns . (5.2) and (5.5) 

0 _ P 0 ^ 1+e) - P 0 
P 1 Lu 2 


so that 


P 2 


1 . . _2 T ,rc 

i. p a. h 


M 2 C 


5.7 


Por a given airplane geometry, € is a function of 
Mach number and angle of attack. The perturbation 6 as a 
function of M and a, should vary least for the best location 
of the static pressure sensing port. This variation is found 
suitable for the nose boon location of the pitot-static tube 
fig. (55-A) . The shape of the pitot port is so designed that 
it ca n sense the pitot-pressure accurately upto - 15° of 


angle of atts.de. 


5.2 Pitot-head, .design f pr_np.se, boon, locations 

Por subsonic flow, the fuselage causes significant 
compression of air ahead of it. This results in a static 
pressure error. The static pressure field ahead of the 
fuselage axis is plotted in fig*(55~B). The pressure 
coefficient falls in a non-linear manner ahead of the body 

In the plane of the symmetry of the aircraft, the 
lateral perturbation velocity is zero. The vertical perturbation 
velocity is created by the flow field on the wing. The 
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expression for 0^ in the plane of symmetry is given by 

- f'Du + \ ((1-M 2 ) u 2 + w 2 ) "i 

0 = t 

P IP # ' 

\ 


.5.8 


Now we go for the design of the contour of the pitot- 

static head to comtensate for the positive static pressure 

perturbation at the static pressure sensing port* To avoid 

abrupt compression, a sufficient length of pitot head is 

provided. Smooth curvature ahead of the static port is given to 

avoid flow s operation. The static port location is taken at 

0.55D ahead of the body where the. static pressure coefficient' 

is found to be 0.1217 for a 45° swept tapered wing-elliptic 

body combination. The body is assumed to be of fineness 

ratio 1:10. At' a distance of 0.55 D forward of the body, 

the value of the wing-hody perturbation static pressure 

coefficient is neutralised by the expansion of the subsonic 

flow over the pitot head at the static port. To make the 

static pressure sensing port insensitive to the angle of attack, 

we choose two static port locations, which are diagonally 

opposite to each other and lie in the lateral direction. 

fig. (55-B ), These static ports are not sensitive to the 

flow direction since they lie. parallel to the flow direction 

created by the wing. The use of a single set of static ports in 
22 

Rosemount tube is a contradicting statement to the use of 

\ 21 

four static ports placed diagionally opposite. The advantage 
of a single set of ports for both subsonic and supersonic flows 
is stated to obviate the need for computation of correction 
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22 

within the central air feta computer. The single set of ports 

22 

chosen for the Rosenount tube are very sensitive to angle of 

21 

attack. The two sets of ports used in the Rosemount tube 
are placed diagonally opposite to compensate for pitching, 
yawing or rolling motions. In our present pitot-static 
head design the use of a single set of diagonally opposite 
static ports is considered to nullify the effects of pitching, 
yawing or rolling motions. The position of the static ports 
is set diagonally opposite in the lateral direction to avoid 
the necessity of using computed corrections within the central 
air data computer. To make the static pressure sensing port 
insensitive to the vertical perturbation velocity and angle of 
attack, we set the value of the vertical perturbation 
velocity created by the pitot-static tube at the static port 
location equal to zero. The mathematical formulation for the 
design of the aerodynamic compensator is described below. 


5 • 3 M athem atical for mulation? 

The effect of pitot - static tube shape is simulated 
by a linearly varying line source singularity. The strength 
of this singularity is determined to get the required shape 
of the pitot-static tube. The potential function due to a line 
source of constant strength located along the x-akis between 
x = o and x = 1 is given by the following expression: 






= m Si 


- 5*3 


r 


(x-O 2 + (1-M 2 )r 2 


5.9 



30 * 


m = strength of source per unit length. She integration is 
carried out between the limits = o and \= 1 and resulting 


potential may* be written as 


x 1 + 


~ Re ^i’" < V' ot i l0 ^ " x“" ay ^ 


...5.10 


whore 


a i..“ + (1-M 2 ) r 2 x 1 = x - x^. 

" 2 = TTl-M^T 2 * 2 - x - 1 


In subsonic flow the distances d^ and d 2 may be 
taken from the origin and tail of the line source to the 

4 

point P (x,rlf’ 1 ) in the Prandtl-Gla$ert coordinate 

system, fig. ( 1-3. ), in supersonic flow d^ and d^ are 

both imaginary ahead of Mach cone from origin; d^ is real 
and d 2 is imaginary between the Mach cone ahead of the body 
and the Mach cone aft of body. Behind the Mach cone both 
d^ and d 2 and are real, the variation of d^ and d^ is shown 
in fig. (1-d). In supersonic flow the distance d^ between the 
Mach cone ahe?,d and aft of body is interpreted as the distance 
from the origin to the point P on the circle of diameter X, • 
fig. (1-d). , 

Ihe logarithmic term in eqn.(5.10) can be expressed, 
as a function of hyperbolic sine and hyperbolic cosine for 
subsonic and supersonic domains respectively. The velocity 
components are given by the partial derivatives of potential 
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function as follows: 

8(j) 

0X 


XL 


V = 


8(j) 

iy 


V/ 


0<j) 


After carrying out the partial differentiation the 
perturbation velocity components can be expressed as 


r x 2 + h . 

u = m Re | log( ) + 

i_ 1 + ~ U 1 


d 2 J 


hug 


n 


- r Re ( 


a i - 
do 


. 5.11 


) = m.d.- 


The pressure coefficient through the Bernoulli l s 
equation for compressible steady flow is given by 

-fun + I C( 1-M 2 ) + V rD 

n . t 5.12 

P c \ U 2 

The desired shape of the nose contour is simulated 
by a pair of line source and line sink singularities 
distributed over lengths 1^ and 1 ^ respectively as shown in 
fig. (55—B) J?or a known length 1 = 1^+ l^j of the nose 
contour, 1^ and 1^ are assumed equal for the sake of simplicity 
and the strengths m^ and m^ of the source-sink singularities 
is determined by eqns. (5.1l) and (5.12). 

The shape of the compensating nose contour which 
provides a negative pressure coefficient equal in magnitude 



to the positive pressure coefficient due to the assumed wing 
"body combination have been calculated and plotted in 
fig. (55-B) . 


5.4 Pitot_head_ design^ for^ fuselage^ dorsal^ mount ingj 

The total perturbation pressure coefficient due to 
the wing~body combination is plotted on the fuselage dorsal 
side in fig. (56-B). The pj_tot static tube is assumed to be 
located at a distance of 0.5P above the top surface of the 
fuselage at the center. She variation of pressure coefficient 
at the static port location is plotted as a function of Mach 
number fig. (56-A) . The pressure coefficient on the fuselage 
dorsal side has a negative value. The shape of the nose 
contour of the pitot-static tube is simulated by a combination 
of two line sources of strengths m^ and m^ and a line sink 
of strength n 2 as shown in fig. (5‘7) The line source 
singularities are distributed over lengths 1^ and 1^, and the 
line sink singularity is distributed over length 1 2 . The 
lengths ±2 and 1^ are expressed as 


1 , 


2 1-j = 21 3 


The sink strength m^ is chosen as half of the 
source strength m^ . The strengths m^ and m^ are determined 
through eqns. (5.11) and (5.12). The equations can be 
expressed as 
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where O n and d^, (rt=1 , 2, 3 ) are defined through eqns.(5.11) 
and (5.12) fig. ( t-i.). The third singularity is a line 
source singularity. Thi_ developes a curved surface aft of 
the static port. This causes a significant compression of 
air on the static port and develops a positive press ur e 
coefficient. To get a smooth shape of the nose contour a 
proper choice of the length over which each singularity is 
distributed, has to be made. Oare has been taken to keep the 
contour curvatures as smooth as possible to avoid abrupt 
compression and expansion of flow past it. The design of the 
pitot-static tube for fuselage dorsal side is complicated. 

The value of the pressure coefficient required to be 

neutralised on the fuselage dorsal side is high and negative 

# 

(-0.1819) as compared to the value of pressure coefficient 
at the nose boom location (+0.1217). This is because the 
pitot-static tube and the wing are centerly mounted on the 
fuselage. The relative location of the pitot-static tube and 
the wing, on fuselage shall be decided to avoid the complicated 
design of the pitot head. The design of the nose contour 
follows a simple and smooth curvature for nose boom location. 

As in the case of the nose boom, the static pressure sensing 
port is made diagnally opposite in the lateral direction for 
the dorsally mounted pitot-static tube. This makes the 
static port insensitive to angle of attack. The nose contour 
of the dorsally mounted pitot-static tube can be computed 
for a given pressure coefficient due to wing-body combination 
by using the programme given in appendix B, 
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Though simple in design, the lengthy pitot-static tube is 
not ^referable, The design of a medium size pitot— static tub© 
is made to provide the negative pressure coefficient to 
neutralise the pressure coefficient existing due to the wing— body 
combination. One disadvantage of this location is apparent 
that the nearer we approach the wing tip, the greater is the 
variation of pressure coefficient with Mach number and , nglc of 
attach. The effect of angle of attack can be readly eliminated 
by the proper location of the static port on the nose contour 
fig. (59-B) . The variation of pressure coefficient with 
Mach number for this location can be neglected over the 
increased length of pitot tube for any other better location 
ahead of the wing tip. The pitot tube designed for wing tip 
location is free from position error for the Mach number 
ranging between (0.3“*0.8). Giving the input data, that is the 
distance of the static port location from the wing tip and the 
value of the perturbed pressure coefficient at this location, 
the design of the tube can be computed by using the programme 
given in appendix B. 



CHAPTER - 6 


6.1 Results^ and^ Discussion; 

The application -of the above method of analysis has 
been done for the steady subsonic flow. The aerodynamic 
compensation extends over "a large range of Mach number. 

Static pressure sensing ports are made insensitive to angle of 
attack. The error in static pressure sensing due to a change 
in angle of attack is least at the nose boom location. The 
static pressure sensing port lies parallel to the plane of 
symmetry and is independent of the flow direction at any angle of 
attack. At very high angles of attack the flow separates 
on the bottom and the topside of pitot-static tube. The 
location of the static port is thus free from seperation of 
flow at any angle of attack. The best location for an 
aero dynamically compensated pitot-static tube is the nose 
boom since the pressure distribution ahead of the nose boom 
does not vary much with Mach number and angle of attack. 

Eor short nosed airplanes to minimise the effect of angle of 
attack, the length of the pitot-static tube should be 
increased since as we go forward of the body the variation in 
pressure coefficient with Mach number and angle of attack 
decreases. 

The variation of pressure coefficient with Mach 
n umb er is more at the fuselage dorsal side as compared to it’s 
variation at the nose boon location. Whore the pressure 
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coefficient lias a high, negative value, a large ' compression of 
air is required to provide the positive pressure coefficient 
at the static port location, for this reason the shape of 
the nose contour consists of double curvature. To avoid 
abrupt compression and expansion the nose contour is made 
smooth. The increased length of the nose contour does not 
effect the overall length of the pitot-static tube. The 
effect of angle of attack is minimised by lateral location of 
static ports. 

In the wing tip location of pitot -static tube the 
pressure coefficient varies significantly with Mach number and 
angle of attack. The pressure coefficient ahead of the 
wing changes steeply with distance. At one location pressure 
coefficient becomes zero. The design of the pitot-static 
tube involves no aerodynamic compensation fo± this location. 

The point of zero pressure coefficient lies at a sufficiently 
great distance from the wing tip. for this reason the length 
of pitot-tube has to be large. This has -disadvantages from. - 
standpoint of maintenance and service, increased weight and 
cost. In addition to a long pitot-static tube an aero dynamic ally 
compensated pitot static tube is designed for a moderate 
length. In the vicinity of the wing tip, the pressure 
coefficient varies significantly with Mach number and angle of' 
attack. This variation decreases with increasing distance 
from the wing tip. The pitot tube for the wing-tip location 
is, therefore, longer than the pitot-static tubes used at the* 
nose boom or fuselage dorsal location. . 
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Aerodynamic compensation and insensitivity to angle 
of attack are the important features of designed pitot-static 
tubes. The mathematical method can- be applied to any airplane 
wing-body configuration, being the computer programming for 
wing-body combination, the value of the pressure coefficient 
can be predicted at any point. With the pressure coefficient 
and the desired length of the nose contour of the pitot-static 
tube as the input data, the computer programme made for Part B 
gives the shapt of the aerodynanic compensator. 
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